Neutron Skin Thickness of 90 Zr Determined By Charge Exchange Reactions 
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Charge exchange spin-dipole (SD) excitations of 90 Zr are studied by the 90 Zr(p, n) and 90 Zr(n,p) 
reactions at 300 MeV. A multipole decomposition technique is employed to obtain the SD strength 
distributions in the cross section spectra. For the first time, a model-independent SD sum rule value 
is obtained: 148 ± 12 fm 2 . The neutron skin thickness of 90 Zr is determined to be 0.07 ± 0.04 fm 
from the SD sum rule value. 
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Proton and neutron distributions are among the most 
fundamental properties of nuclei. Proton distributions 
are precisely known from the charge distributions deter- 
mined by electron scattering 0]. On the other hand, 
our knowledge of neutron distributions, which have been 
studied mainly by hadron-nucleus scattering, is limited 
because descriptions of strong interactions in nuclei are 
highly model-dependent 0. Reliable neutron distribu- 
tions will improve the understanding of the nucleus and 
nuclear matter 0, 0, IE IE ■ Recent theoretical stud- 
ies using the Skyrme Hartree-Fock (HF) and relativistic 
mean-field models 0, IE] have shown that the neutron 
skin thickness, defined as the difference between the root 
mean square (rms) radii of the proton and neutron distri- 
butions, imposes a strict constraint on the neutron mat- 
ter equation of state, which is an important ingredient 
in studies of neutron stars [E ■ It is also known that 
the neutron skin thickness is strongly correlated with the 
nuclear symmetry energy 0- Reliable neutron dis- 
tributions are also needed for analyses of atomic parity 
violation experiments [id ITlT ] and of pionic states in nu- 
clei H2. 

Several attem pts have been made to determine neutron 
distributions 0, EE EE EE EE] • Ray e * a l- analyzed pro- 
ton elastic scattering on several nuclei at 800 MeV using 
impulse approximation and obtained a neutron thickness 
of 0.09 ± 0.07 fm for 90 Zr 13]. The cross sections for 
excitation of isovector giant dipole resonances with al- 
pha scatterings were measured at KVI and neutron skin 
thicknesses of 116,124 Sn and 208 Pb were obtained EE] with 
uncertainties of ±0.12 fm. Trzcihska et al. measured the 
strong-interaction effects on antiprotonic x-rays on vari- 
ous targets ranging from 16 to 238 U and obtained a neu- 
tron skin thickness of 0.09 fm for 90 Zr with a statistical 
error of ±0.02 fm EE]- Unfortunately, these analyses are 
model-dependent. Parity violation electron scattering is 
a promisin g to ol for observing the neutron distributions 
cleanly [IEEE]> although no data are available so far. 

An alternative method for determining the neutron 
rms radius is provided by the model-independent sum 
rule strength of charge exchange spin-dipole (SD) excita- 
tions EH- The operators for SD transitions are defined 



by 



S± = 



with the isospin operators ts = t z , t± = t x 
model-independent sum rule is derived as 

S--S + = ^-(N(rX-Z{r\) 



± ity. 



The 



(2) 



where S± are the total SD strengths. The mean square 
radii of the neutron and proton distributions are denoted 



as 



and (t 



p , respectively. Thus, the rms radius 
of the neutron distribution yj (r 2 ) n or the neutron skin 
thickness S np — \J (r 2 ) n — y/ (r 2 ) p can be derived from 
Eq. J5J by using the rms radius of the proton distribution 



obtained from the charge radius if the sum rule 



value (S- — 5+) is obtained experimentally. 

To obtain the SD sum rule value, Krasznahorkay et 
al. measured the ( 3 He, t) reaction on tin isotopes at 450 
MeV j2ij. The cross section spectra at 6 = 1° were ana- 
lyzed by a peak fitting technique assuming the Lorentzian 
shape, and the S— value was obtained. Since there was no 
(n,p)-type measurement, they determined the S+ value 
by assuming the energy-weighted sum rule in a simple 
model where the unperturbed particle-hole (ph) energies 
are degenerate at a certain energy E3 • The overall nor- 
malization of the sum rule value was done by using the 
calculated result of neutron skin thickness of 120 Sn. To 
determine the sum rule value model-independently, how- 
ever, one needs to perform (n,p)-type experiments as well 
as to identify the AL = 1 cross sections in the observed 
spectra. Then, one needs to relate the AL = 1 cross 
sections to the SD strengths. 

This paper represents the first attempt to extract the 
model-independent sum rule value from both (p, n) and 
(n,p) reactions on 90 Zr. The dipole components of the 
cross section spectra are identified by multipole decom- 
position (MD) analysis [2lj] of the 90 Zr(p, n) and 
90 Zr(n,p) H3] data at 300 MeV. It should be noted that 
at 300 MeV, the spin-flip cross sections are large, while 



2 



the distortion effects are minimal [24|. Thus, the charac- 
teristic shapes of the angular distributions for each an- 
gular momentum transfer (AL) are most distinct. The 
contribution of the non-spin dipole component is small 
due to the energy dependence of the effective interaction. 

The MD analyses were performed on the (p, n) and 
fn, p) excitation energy spectra in a consistent man- 
ner [2I], [23j to obtain various AL components of the cross 
section. We use here the extracted AL = 1 cross section 
spectra given in Fig. 2 of Ref. [2]|. The SD strengths 
are obtained by assuming a proportionality relation sim- 
ilar to that established for Gamow- Teller excitation. The 
proportionality relation between £?(SD) and the AL = 1 
component of the cross section, <tal=i ±(q, oj), is given 
by 



CAL=i,±(g,w) = CT S D ± (g,^)S(SD ± ) 



(3) 



where ctsd ± (<Z,w) is the SD cross section per £?(SD±) 
and depends on the momentum transfer q and the energy 
transfer to. The <tal=i(<Z, w) data for the (p, n) and (n,p) 
channels were taken from the result of MD analysis at 
4.6° and 4-5°, respectively. The <7al=i(<?, w) spectra at 
these angles are most sensitive to the SD cross sections 
since the corresponding momentum transfers are q = 0.3- 
0.4 fm _1 and thus, the AL = 1 cross sections take on 
maximum values. 

Distorted wave impulse approximation (DWIA) calcu- 
lations are used to obtain <tsd ± (?) u )- The DWIA calcu- 
lations are performed with the computer code DW81 |2^| 
for the 0~, 1 — , and 2~ transitions. The one-body tran- 
sition densities are calculated from pure lplh configura- 
tions. All lhw configurations were examined. The op- 
tical model potential (OMP) parameters are taken from 
Ref. [26j . The effective NN interaction is taken from the 
t-matrix parameterization of the free NN interaction by 
Franey and Love at 325 MeV [27J. The radial wave func- 
tions are generated from a Woods-Saxon potential |28j . 
adjusting the depth of the central potential to reproduce 
the binding energies. Details of the calculations are found 
in Refs. plll^. 

The calculated SD_ unit cross sections at 6 = 4.6° 
and w = MeV are 0.28 ± 0.03, 0.24 ± 0.06, and 
0.29 ± 0.05 mb/sr/fm 2 for the 0", 1", and 2" transi- 
tions, respectively. The uncertainty indicated for each 
multipole is the dependence of the cross sections on the 
ph configurations. Gaarde et al. pointed out that the 
tensor term of the effective NN interaction affects the 
proportionality among the multipoles in the study of SD 
excitations in the 12 C(p, n) reaction at 160 MeV. At 300 
MeV, however, the ratio of the amplitude of the isovector 
tensor interaction to that of the isovector spin (err) in- 
teraction at q — 0.3-0.4 fm -1 is smaller |2jj so that the 
calculated unit cross sections are close to one another. 
To check the validity of the above calculations with pure 
lplh configurations, the calculations of SD unit cross sec- 
tions are performed by using the transition densities for 
several states obtained by the HF+RPA (random phase 
approximation) calculations |30j |. The calculated unit 



20 



> 

S 10 



d> 
!h 

to 10 

P 

in 



dfl(SD_) 
dE 


• ' 1 ' ' ■ ' 1 ' ' ' ' | ' ' ' ' 1 ■ 

• 90 1-7 / \ 

~ ♦ Zr(p,n) 

m \, Drozdz et al. 

/ *\ 


/V, 




/ • 








. . 1.^*1* . 


■ . I . . . . I . . . . I . . 


dfl(SD + ) 
dE 


i 90 Zr(n no) 




(shifted by +17MeV) - 

..■ V1T, '\itTll* 



10 



20 30 
E x (MeV) 



40 



50 



FIG 

and 



1: Charge exchange SD strength — — - (upper panel) 



dB(SD + ) 
dE 



(lower panel). The circles and squares are 

the experimental data. The d " B ^? + ^ spectra are shifted by 
+17 MeV. The curve is the results of the second RPA calcu- 
lation by Drozdz et al. [3^|. 



cross sections are found to be consistent with the above 
values. 

Since the main subject of this study is the total SD 
strength, rather than the individual strength of each 
transition, we use the averaged unit cross section of 
(T S d_(4.6°,0 MeV) = 0.27 mb/sr/fm 2 in the analysis. 
Similarly, the unit cross section in the SD + channel is 
estimated to be <7 SD+ (4-5°, MeV) = 0.26 mb/sr/fm 2 . 
The unit cross sections are calculated at each energy bin 
of the cross section histogram. The systematic uncer- 
tainty in <tsd ± (<?, w) due to the input parameters for 
DWIA calculations has been evaluated by using har- 
monic oscillator radial wave functions and by other sets 
of OMP |24i |2|| • The systematic uncertainty thus esti- 
mated is 14%. 

The dB d ^ D ^ distributions obtained by using Eq. © are 
shown in Fig. ^ The horizontal axis in the dB<y ^ ^ 
spectrum is the excitation energy of the residual 90 Nb 
nucleus. The SD_ strength spectrum shows a dominant 
resonance structure centered at E x = 20 MeV and the 
strength extends to ~ 50 MeV excitation. The — Iffi^ 
distribution is shifted by +17 MeV to account for the 
Coulomb displacement energy |28j and the nuclear mass 
difference. The SD + strength distribution forms a broad 
bump centered at 28 MeV (or 11 MeV in 90 Y) with a 
width of - 15 MeV. 

Drozdz et al. studied the SD_ strengths with the RPA 
model including the coupling between lplh and 2p2h 
states, though the calculation is not self-consistent, us- 
ing a phenomenological Woods-Saxon potential to ob- 
tain the ground state wave function |33|. They found 
that the mixing of 2p2h states results in a large asym- 
metric spread in the strength of the SD resonances, with 
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FIG. 2: Integrated charge exchange SD strengths. The upper 
panel shows the S- and S+ spectra. The lower panel shows 
the S- — S+ spectrum. 



about 30% of the total strength shifted to excitation en- 
ergies above 28 MeV j3i|. The curve in Fig. 1 shows the 
predicted strength distribution. This calculation gives 
a reasonable description of the strengths above 20 MeV 
although it overestimates the strengths at lower excita- 
tion energies below 15 MeV. To discuss the details of 
SD strength distributions and their relations to the neu- 
tron skin thickness, the self-consistent HF+RPA models 
should be employed where the same two-body interaction 
is used throughout the calculations. 
The integrated SD strength, 



S± 



dB(SD d 
dE~ 



-dE, 



(4) 



is plotted in Fig. [21 While both integrated SD strengths 
increase steadily, the sum rule value, S*_ — S+, remains 
almost constant in the excitation energy range of 30- 
50 MeV, where it lies within 149 ± 5 fm 2 . Since the 
MD analysis in the (n,p) channel is unstable above an 
excitation energy of 40 MeV in Fig. ^ we integrate the 
strengths up to 40 MeV. The SD_ strengths are 5_ = 
247±4(stat.) ± 12(MD) fm 2 , where the statistical uncer- 
tainty and the uncertainty in the MD analysis are given. 



The corresponding S + value is 98±4(stat.)±5(MD) fm 2 , 
integrated to an excitation energy of 23 MeV in 90 Y 
(40 MeV in Fig. 01. The sum rule value yields S- -S+ = 
148 ± 6(stat.) ± 7(syst.) ± 7(MD) fm 2 , where the system- 
atic uncertainty of the normalization in the cross section 
data (5%) is also included. The 14% uncertainty in the 
SD unit cross section is not included. 

If this sum rule value is interpreted in terms of Eq. , 
we obtain N{r 2 ) n - Z(r 2 ) p = 207 ± 17 fm 2 , where the 
statistical and systematic uncertainties are combined in 
quadrature with the uncertainty in MD analysis. The 
rms radius of the proton matter in 90 Zr is estimated to 
be 4.19 fm after correcting for the effect of the proton 
form factor from the charge radius |32| . The neutron skin 
thickness and the rms radius of the neutron distribution 
calculated from Eq. (J2J are S np — 0.07 ± 0.04 fm and 
\J (r 2 ) n — 4.26 ± 0.04 fm, respectively. The S np value 
obtained in the present study is consistent with that 
obtained from the analysis of proton elastic scattering 
(0.09 ± 0.07 fm) [13 but with a smaller uncertainty. Our 
results also agree with the neutron skin thickness deter- 
mined by antiprotonic x-ray measurements by Trzcihska 
et al. (0.09 ±0.02 fm) [j^. though their analysis contains 
some assumptions whose uncertainties are not specified. 
We note that the accuracy of \J {r 2 ) n in this study is 1% 
level, which is the same as the goal of the parity viola- 
tion experiment of electron scattering at Jefferson Lab- 
oratory [18j. To improve the reliability of the present 
analysis, the SD unit cross sections should be studied in 
the mass region around 90 by measuring the cross sec- 
tions of the SD transitions whose strengths are known 
from decay measurements. 

In summary, a consistent MD analysis of the (p, n) and 
(n, p) reaction data from 90 Zr has been performed and SD 
strength distributions of both channels are obtained ex- 
perimentally for the first time. The integrated strengths 
obtained experimentally are S- — 247 ± 4(stat.) ± 
4(syst.)±12(MD) fm 2 and S+ = 98±4(stat.)±4(syst.)± 
5(MD) fm 2 , up to 40 MeV and 23 MeV, respectively. By 
using the two experimental sum rule values, the model- 
independent formula (J2J yields a N(r 2 ) n — Z(r 2 ) p value 
of 207 ±17 fm 2 , which corresponds to a neutron skin 
thickness of S np = 0.07 ± 0.04 fm. The method used in 
this work is applicable to heavy and medium heavy nuclei 
on which both the (p, n) and the (n,p) measurements are 
possible. 
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